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A.O  EQ SOLVE  Program  Documentation 
A.l  Introduction  — Program  Overview 
A. 2 Detailed  Main  Program  Structure 
A.  3 Program  Calculations 
A. 4 Detailed  subroutine  structure  — GRIM 


A. 0 EQSOLVE  PROGRAM  DOCUMENTATION 

This  appendix  describes  a software  program  called  EQSOLVE,  which 
was  developed  at  TRW  Defense  and  Space  Systems  Group  for  the  ONR  in  con- 
junction with  this  contract.  EQSOLVE  accepts  parameters  of  the  Choma  Bipolar 
Transistor  model  and  a tabular  list  of  base  current  values,  and  produces 
device  performance  curves  such  as  Beta  vs.  Ic.  The  performance  curves  are 
presented  in  families,  allowing  parametric  variation  with  a third  parameter 
which  may  be  a model  parameter  or  the  collector-emitter  voltage.  This 
allows  a sensitivity  analysis  of  transistor  operating  characteristics  with 
respect  to  model  parameters  and  allows  a "fine  tuning"  of  the  model  based 
on  easily  observed  device  performance  curves.  The  EQSOLVE  code  is  made 
to  accept  model  data  from  a companion  code  called  PAREV.  PAREV  generates 
from  electrical  and  processing  data  the  model  parameters  used  as  input 
to  EQSOLVE. 
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A. 1 INTRODUCTION  - PROGRAM  OVERVIEW 

A. 1.1  Objectives 

The  EQ SOLVE  program  is  a dual-purpose  code.  In  the  early  stages 
of  contract  performance,  the  code  provided  insight  into  development  of  test 
methodology  for  specific  model  parameters.  By  examination  of  Calcomp  plot 
families  of,  for  example.  Beta  vs.  Ic  with  extrinsic  emitter  resistance  Rg 
as  a parameter,  the  effect  of  model  parameters  on  test-generated  curves  is 
made  evident,  and  insight  is  gained  into  testing  methods  for  various  model 
parameters.  In  the  later  stages  of  contract  performance,  after  parameter 
testing  methodology  has  been  established,  the  EQSOLVE  is  code  is  used  to 
fine  tune  and  to  verify  the  performance  of  a given  device  model  prior  to 
model  release  for  CAD  activity.  An  additional  future  goal  for  the  EQSOLVE 
is  in  the  area  of  computer-assisted  device  design,  i.e.,  the  process  of 
using  an  ideal  set  of  performance  characteristics  to  direct  the  device 
fabrication.  EQSOLVE  provides  a link  allowing  software  comparison  of  ideal 
and  predicted  device  performance  in  an  iterative  improvement  feedback  loop 
that  can  produce  masking  and  diffusion  information  for  device  fabrication 
based  upon  desired  performance  characteristics. 

A. 1.2  Program  Structure 

The  EQSOLVE  program  is  written  in  FORTRAN  IV  for  use  on  the 
CDC  6000  or  CYBER  series  machines.  The  core  requirements  are  less  than 
60K  (octal)  to  insure  interactivity  capability.  Expected  interactive 
terminal  time  for  a single  graph  output  is  less  than  five  minutes  with  a 
total  cost  of  less  than  $5. 

A multiple  graph  output  capability  exists,  allowing  a user  to 
generate  several  performance  curve  families  on  which  a different  model 
parameter  is  varied  for  every  family.  Additionally,  different  performance 
curves  for  a given  model  parameter  (e.g, , f vs.  I with  R varying  and 
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B vs.  I£  with  R£  varying)  can  be  generated.  Each  family  request  generates 
a list  of  coamon  device  values  in  addition  to  the  Calcomp  plot,  namely: 


t 


i 


1.  Ig  — Base  current 

2.  Ic  — Collector  current 

3.  hpE  — DC  current  gain 

4.  Bp  — Base  pushout  factor 

5.  fT  — Gain  bandwidth  product 

6.  VgE  — Base-emitter  voltage 

Program  input-output  interaction  and  user  interaction  are  sum- 
marized in  figure  A-l.  User  interaction  consists  of  real  time  requests, 
wherein  the  details  and  contents  of  the  Calcomp  families  of  curves  are 
stipulated. 


TAPE  2 

l 

TAPE  4 


TAPE  5 


Cal  comp 
family 


TAPE  7 


Figure  A-l  EQSOLVE  Input-Output  Interaction 
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DETAILED  MAIN  PROGRAM  STRUCTURE 
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A. 2.1  Namelist  Input 

The  namelist  input  works  in  a default  override  fashion.  Default 
parameters  corresponding  to  the  model  of  a 2 Gigahertz  small  geometry 
bipolar  transistor  are  used  on  input  TAPE5.  The  program  reads  the  default 
values  and  then  reads  in  any  user-supplied  update  values  on  input  TAPE2. 
The  TAPE2  values  replace  the  default  values  for  the  duration  of  the  user 
session.  If  the  PAREV  companion  code  is  used  prior  to  EQSOLVE,  TAPE2  is 
available  with  a complete  set  of  updated  parameters,  and  it  totally  over- 
writes the  TAPE5  default  set.  This  is  the  expected  mode  of  usage.  An 
example  TAPE2  for  PAREV  output  and  EQSOLVE  input  is  shown  in  figure  A-2. 
Detailed  parameter  explanation  is  found  in  Section  2 of  this  report,  and  a 
list  of  input  parameters  with  default  values  is  supplied  in  table  A-l. 

A. 2. 2 Base  Current  File  Input 

Base  current  file  TAPE4  contains  the  list  of  base  current  points 
at  which  the  bipolar  device  performance  is  to  be  calculated.  The  user 
thus  has  flexibility  in  selecting  the  range  and  coarseness  of  Calcomp 
output  by  tailoring  TAPE4  to  his  own  specifications.  Input  format  is 
free-format  as  follows: 

N < integer  value  “ number  of  base  current  points  > 

IB1  First  base  current  point  in  amperes 

IB2  Second  base  current  point  (1^  > 


!l 


! 


IBN  Nth  base  current  point 

An  example  of  a TAPE4  file  is  given  in  Figure  A-3.  Base  current  points 
are  from  .luA  to  1.0mA. 
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C LIST,TAPE2 
1ITRANS 


! 

RBO 

3 

1.8, 

RBB 

3 

1.3E+02, 

RE 

= 

5.0E-01, 

RCC 

3 

0. , 

* 

BS 

= 

3.333333333E-01, 

BC 

= 

3 • 333333333E— 0 1 , 

\. 

BE 

= 

5.0E-01, 

f I 

CSO 

3 

2.0E-12, 

* 

CEO 

= 

8.998858595E-13, 

CCO 

7.999348551E-13, 

PHISS 

= 

9.0E-01, 

PHICO 

3 

9.004004474E-01, 

PHIEO 

3 

1.200517362, 

ISR 

= 

3.0E-12, 

15 

= 

3.2E-16, 

IER 

3 

3.0E-14, 

! 

IBR 

= 

3.0E-13, 

VER 

= 

2 . 1 96833685E+01 , 

' 

VCR 

3 

6.056999515E+01, 

1 

NC 

3 

1.8, 

VT 

3 

2.587500647E-02, 

VTS 

= 

3.0E-02 , 

OFN 

3 

6.178E-16, 

8RN 

= 

2. 26141181 2E-12, 

TAUFO 

3 

1.273024046E-10, 

TAUR 

= 

4.667132469E-07, 

BETAO 

= 

5.5E+01, 

BE  TAR 

3 

1.5, 

AE 

3 

9.8E-07, 

1 

XO 

3 

2.7E-04, 

UCPRIHE 

3 

3.0E-04, 

NE 

3 

1.5, 

KR 

3 

5.0E-01, 

# 

BF2 

3 

1.4E+01, 

VB 

3 

1 .3681 20028E+01 , 

VF 

3 

2.587500647E-02, 

FF 

3 

I.027800526E-14, 

VCE 

3 

5.0, 

IBB 

3 

1 .380000343E-03, 

PC 

3 

3.25E-03, 

BEND 

Figure  A-2  PAREV  Output  - EQSOLVE  Input  TAPE2 
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TABLE  A-l  INPUT  PARAMETER  LIST  - SEGMENT  1 OF  3' 


RCC  RCC  Low  Injection  Extrinsic  7.0  ohms 

Collector  Response 


TABLE  A-l  INPUT  PARAMETER  LIST  - SEGMENT  2 OF 


Nonideal  Emitter  Exponent 
Multiplier 


TABLE  A-l  INPUT  PARAMETER  LIST  - SEGMENT  3 OF  3 
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• IE  — 6 

• 2E-b 

• 4 E—6 
.7E-6  ' 
1.0E-6 
2.0E-6 
4,0c-b 
7 .OE-b 
10.E-6 
20.E-6 
40.E-6 
70.E-6 

• 1 E-3 

• 2E-3 

• 4E  — J 

• 7E-  3 
i.e-3 


Figure  A-3  TAPE4  - Base  Current  Input  File  Example 


A. 2. 3 Tabular  Output 

Local  file  TAPE7  is  written  at  every  Ig  point  with  the  variables 

Vc>  Ig,  Ic>  hpg.  Bp,  f,p,  and  Vgg.  TAPE7  may  be  listed  at  the  termination 
of  an  interactive  session  to  supplement  the  Calcomp  printer  plots . An 
example  TAPE7  partial  listing  is  shown  in  Figure  A-4. 

A. 2. 4 Calcomp  Output 

Examples  of  Calcomp  output  graphs  are  shown  in  the  sensitivity 
curve  sets  of  Appendix  B. 

A. 2. 5 User  Interaction  - Sample  Session 

The  interactive  session  using  EQSOLVE  requires  user  input  in 
response  to  terminal  prompters  for  specification  of  Calcomp  output 
details.  The  first  request  identifies  the  independent  (x-axis),  dependent 
(y-axis) , and  parametric  variables  by  name  and  index  number.  Index 
number  is  the  relative  position  of  the  variables  within  a COMMON  block 
called  PRAMTRS,  and  is  used  for  Internal  program  identification  only. 
Parameter  name  is  used  for  axis  labeling  and  may  be  specified  by  alpha- 
numeric string  less  than  eight  characters  without  altering  graphical 
output.  Index  values  are  tabulated,  including  default  values,  in 
Figure  A- 5. 
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C LIST- TAPE? 


VCPRIME  * -4.572E+00 

IB  = 1 . 000E-07 

I C = 1 .£76E— 06 

HFE  = 1 .276E+01 

VEE  = 5.713E-01 

EF 

= l.OOOE+GO 

FT 

— 4.639E+06 

VCPRIME  = -4.554E+00 

IE  = £ . 000E-07 

IC  = 3.£l£E-06 

HFE  = 1 .606E+01 

VEE  = 5.953E-01 

BF 

= 1 . 001E+0Q 

FT 

= 1 . 143E+07 

VCPRIME  = -4.536E+G0 

IB  = 4 . 000E-07 

IC  = 7.390E-06 

HFE  = 1 .973E+01 

VEE  = 6.136E-01 

EF 

= 1.003E+0G 

FT 

= £ .719E+07 

VCPRIME  = -4.5££E+00 

IB  = 7 . G00E-07 

IC  = 1.600E-05 

HFE  = £.£86£+Gl 

VEE  = 6.37GE— 01 

EF 

= 1.006E+0G 

FT 

= 5.298E+07 

VCPRIME  = -4.514E+0G 

IB  = 1.000E-06 

IC  = £ .491E-05 

HFE  = £ .491E+G1 

EF 

= 1 . 009E+00 

FT 

= 7.939E+G7 

VEE  = 6.485E-01 


Figure  A-4  TAPE7  Partial  Output  List 
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X>LOK>hJ0'OCX)^IONLn-t'LOKJ»-»OvOOO'^lONLn^U)K>h-*Ov£>OO^JONUi4>U>ro»-* 


Parameter 

Default 

Index 

Parameter 

Default 

RBO 

9. 

35 

NDC 

1.E16 

RBB 

150 

36 

PHIC 

.7 

RE 

.5 

37 

DNC 

48. 

RCO 

5. 

38 

XO 

1.2E-4 

PC 

.00325 

39 

KR 

.4 

BS 

.33333 

40 

TD 

20.E-12 

BC 

.33333 

41 

M 

2.0 

BE 

.5 

42 

NE 

1.5 

cso 

2E-12 

43 

IB 

CEO 

•9E-12 

44 

VCPRIME 

CCO 

.8E-12 

45 

IBB 

PHISS 

.9 

46 

QF 

PHICO 

.9 

47 

QR 

PHIEO 

1.2 

48 

LAMBDAE 

ISR 

3.E-12 

49 

ISS 

IS 

3.2E-16 

50 

IN 

IER 

3.2E-14 

51 

II 

IBR 

3E-13 

52 

10 

VER 

20. 

53 

ILIM 

VCR 

60. 

54 

BF 

NC 

1.8 

55 

QBN 

VT 

.0259 

56 

VBE 

VTS 

.03 

57 

IC 

QFN 

1.5E-16 

58 

HFE 

QRN 

1.5E-12 

59 

FT 

TADFO 

159.E-12 

60 

VCE 

5.0 

TAUR 

1270E-12 

61 

BETAE 

BETAO 

55. 

62 

BF2 

13.64 

BEXAR 

1.5 

63 

VB 

17.793 

W 

.6E-4 

64 

VF 

2.59E-3 

AE 

9.8E-7 

65 

FF 

1.336E-14 

VLIM 

6.E6 

66 

F 

DNB 

22.6 

67 

BFO 

WCPRIME 

3.E-4 

Figure  A-5  Common  Parameter  Indices 


A typical  prompter-response  dialogue  is  as  follows: 

ENTER, IDNKE, INDEX, DPNME, INDEX , PARNME , INDEX 

? |*VBE*,56,*IC*,57,*RBO*,l  fr-  user  supplied 

The  dialogue  requested  an  IC  versus  VBE  plot  with  resistor  RBO  as  the 
varying  parameter. 

The  next  dialogue  requests  Calcomp  type,  number  of  curves  per 
plot,  and  the  varying  parameter  values.  Curve  type  is  specified  by  the 
following  index  list: 


X-Axis 

Y-Axis 

Index 

LINEAR 

LINEAR 

1 

LINEAR 

LOG 

2 

LOG 

LINEAR 

3 

LOG 

LOG 

4 

A prompter- response  dialogue  requesting  a linear  (x)  - log  (Y)  plot  with 
three  values  of  a varying  parameter  equal  to  90,  9,  and  .9  (RBO  values) 
is:  (NCURVES  < 5) 

ENTER, TYPE , NCTJRVES , PARAM  VALUES 
?|2,3,90.,9., .9H-  user  supplied 

The  third  and  final  dialogue  requests  adjustments  in  size  of  the  x-axis 
and  y-axls,  and  number  of  points  between  sumbols  on  a curve.  A sample 
dialogue  requesting  a 6 inch  x-axis , 7 inch  y-axis,  and  plot  symbols  every 
3 points  (default  values)  is: 

UPDATE  XINCH , YINCH , MARKS ( 5 ) IF  REQD 
SCURV 

| YINCH-7  ,X INCH-6  .MARKS" 3, 3 . 3$  user  supplied 

Note  that  three  entries  of  three  (3)  are  required  since  three  curves 
are  requested  by  NCURVES  prompter. 

An  actual  teletype  session  for  Calcomp  output  using  the  above 
dialogues  is  given  in  Figure  A-6.  The  def atilt  x-axis , y-axls  size  and 
symbol  density  are  used  by  supplying  a terminating  dollar  sign  ($) . 

An  actual  Calcomp  output  for  this  session  is  given  in  Figure  A-7. 
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07/13/77.  10.37.10. 

TRU/TS8  Cl  74  - 07/11/77. 

USER  NUMBER i 
JOB  NAME  PLR6. 

PORT  220. 

C GET ,TAPE4,EQS0LVE,TAPE5sEGDATAi 
t RX,I*EQSOLVE  — 

C LGO 

ENTER  INDNNE , INDEX, DPNHE , INDEX, PARNNE , INDEX 
? *VBE*,54t*IC*,57,*RB0*f  1 
ENTER  TYPE,NCURVES,PARAM  VALUES 
T 2 3 90  9 9 

UPDATE  YINCH,XINCH,HARKS<5)  IF  REQD 
9CURV 
? « 

TYPE  1 FOR  ANOTHER  CASE 
T 1 

ENTER  INDNNE , INDEX ,BPNHE , INDEX, PARNNE , INDEX 

<etc.> 


Figure  A-6  Interactive  Session  Printout 


rr 


A. 3 


Program  Calculations 

The  EQSOLVE  program  performs  calculations  of  Internal  and  external 


parameters  of  a given  modeled  transistor  at  a set  of  base  current  values 
and  at  a fixed  VCE  in  the  normal  (forward)  regime.  Due  to  the  dependencies 
of  the  model  equations,  a nonlinear  root  finder  routine  is  required  to 
calculate  internal  currents  and  voltages,  and  an  Iteration  path  is  used 
to  converge  the  internal  currents  and  junction  voltages  to  values  which 
satisfy  both  the  model  equations  and  the  external  constraints  of  a fixed 
base  current  IB  and  collector-emitter  voltage  VCE.  An  overview  of  the  program 
flow  is  shown  in  Figure  A-8,  illustrating  the  internal  variables  Z,  V „ 
and  Ig  generated  by  three  separate  calls  to  a root  finder  routine  called 
GRIM.  An  auxiliary  subroutine  is  used  in  each  call,  which  describes  the 
controlling  equation  containing  the  variable  to  be  found.  Details  of  the 
GRIM  and  auxiliary  subroutines  are  explained  in  other  sections  of  this 
document. 

The  set  of  calculations  begins  with  a requested  current  point, 

Ig,  and  model  parameters  (user  supplied  or  default) . An  Initial  guess  of 
collector  junction  intrinsic  voltage  is  given  as 


V'  = -V 
C CE 


and  GRTM  is  called  to  evaluate  Z from  the  equation 


(A-l) 


Given  Z, 


ZTan(Z)  - yy 


R3B[Tan(Z)  - Z] 
4*Z*Tan(Z) 


(A- 2) 


(A- 3) 


and 


Agp  - AE*Sin(2Z)/2Z 


(A-4) 


The  estimate  V'  is  now  required  to  compute  V'  from 

C £ 


L 


A-16 


h ■ V“p<te/neV  - 1]  + — fe*p<VV  - 13 

o 


+ i^«P<W  - 1]  + I8R[exp(V'/NcVT)  - l] 


(A-5) 


by  another  call  to  the  root  finder  routine.  The  intrinsic  junction 
voltages  are  now  known,  and  the  current  and  charge  terms  are  calculated. 
Prior  to  this,  if  is  not  given  as  input,  the  calculation 


is  performed.  Then, 


(A-6) 


In  - V«P<W  - 1] 

(A-7) 

Ij  - Is[exp(V'/VT)  - 1] 

(A-8) 

Qp  " Qpn  [exp(V'/VT)  - 1] 

(A-9) 

\ ■ «„[«pcw  - 1] 

(A- 10) 

XSS  “ ISR  [«Kp(Vc/VTSJ  - l]  . 

(A- 11) 

calculations  is  Early  effect  and  base  pushout 

parameters . 

XB  • 1 + V'eR  + Wa.  • 

(A-12) 

F - FF  [exp  - l]  , 

(A-13) 

BF0  - 1 + * 

(A-14) 

®F  " ®F0  + (1  " BFO)C(exp(‘F)J  • 

(A-15) 

A-18 


The  third  call  to  the  root  finder  routine  is  made  to  evaluate  Ic  from 
the  equations 


qb»i  ■ r + [(r)2  + bfqf  + qr]1/2 

(A-16) 

xc  “ 

Is 

ISS  - IBRtexP(V6/NCVT)  - 1]  - e;  - 1] 

+ \'h 
qbn 

(A-17) 

Additionally,  the  program  evaluates 

. *o/BF-1\1/2 

CC  V »F2  / 

(A-18) 

Rc  " Rcc(1  " xcc) 

(A- 19) 

'be  • r c«p<w  - 11 

o 

(A-20) 

Ibc  “ Cexp(Vc/VT)  - l] 

(A-21) 

VBE  ' V*BO  + (IBE  + hc>'h  + Vi+  (IB  + V’h 

(A-22) 

Finally, 

'CE  ' VX  ■ VRCO  + «C  - VRC  - V6  + V£  + (IB  + V'h 

(A- 2 3) 

The  comparison  of  to  is  now  performed. 

a 

If  |vc£  - V^l  <_  .05  V^,  the  calculations  proceed.  Otherwise, 
a new  estimate  of  V'  is  generated  by 


and  program  contr^.  s transferred  to  the  second  call  to  GRTM  following 
equation  (A-5) . 


Once  a "suitably  close"  value  for  V__  is  found,  the  program 
calculates  several  parameters  of  interest  at  the  operating  point  (IR, 

V ).  The  first  quantities  are  the  base  pushout  and  forward  base  charge 
sensitivities, 


(Bfq  - l)-FF-exp(V£/VF  - F) 
SBE  " Bf*Vf 

BF2-exp(V'/VB)«[l  - exp(-  F) ] 

SBC  ‘ 

[IC,Bf‘Qfn/(Is‘Vt)1(1  + VT'SBE) 
* " 2*QBN  - XE 


(A-25) 


(A-26) 

(A-27) 


SQG  ‘ (2‘Vcr‘Qbn  +2*Vcr/(2.Qbn  - V^BN  + 2*VCR*  Wc^Bc/V  (A‘ 


The  forward  and  reverse  incremental  transconductances  are  given  by 


gmf  - (IC/VT)-(1  - SQ£/VT)  (A-29) 

gmt  - IC.SQ(.  (A- 30) 

The  program  now  checks  if  base  pushout  is  in  evidence  (Bj,  > 1) 
and  calculates 


Roc- 


Rcc 


Rc- 


C KCC  ^ gmf  gmr  J 

2*Bf2*[(Bw  - 1)/\?]1/2 


(A-31) 


bf>  1 


F2- 


and 


R, 


OB 


2*vt/(3*iBb) 


(A- 32) 


28) 


A-20 


The  junction  capacitances  are  assigned  values  by  the  equations 


(i  - VVbe 


VVF 

cro*a  + 4-1) 


*E0 


v:.  > * 

E ?E0 


(A-33) 


w - VW  c 

Bc'vc 

cco-<1  + -r7> 

^CO 


Diffusion  capacitances  are 


VC  < ° 


V'  > 0 


(A- 34) 


CDE  " + W 


(A- 35) 


CDC  * Bf’TFO*  ^C^BC  ” 


(A-36) 


yielding  total  capacitances  given  by 


CTE  “ CDE  + CjE 


(A-37) 


CTC  “ CDC  + CjC 


(A- 38) 


Substrate  voltage  is  assumed  equal  to  collector  junction  voltage,  and 
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ss 


'so 


(l  - v-/*ss)  s 
cso*(1  + WW 


vc  < 0 


vc  > 0 


(A- 39) 


The  time  moment  calculations  proceed  with 


*BE  1 


8be  " vT  + nE-vT  *[exp(vE/nE-vT)  - l] 


(A-40) 


Fe  - 1 + (gmf  + g^).^ 


8b  * 8be^fe 


CA-41) 

(A-42) 


CE  " CTE/,FE 

(A-43) 

go  - gmr/F£ 

(A- 44) 

gm  - (gmf  - gmr)/FE 

(A-45) 

Roc  + Rco 

\ * l+go.(R0C+RC0) 

(A-46) 

rll  “ 1/gB 

(A-47) 

r22  " K + (1  + gm*RL)*rll 

(A-48) 

r33  “ R0B  + rll  + (1  + 8m*rll),R{ 
<1  + 8°,R0C) *RC0 

^ " 1 + go*(R0C  + Rco> 


and 


T1  " rll*CTE  + ~ Kr)  *r33*cxc  + r44’CSS, 


(A-49) 

(A-50) 

(A-51) 


A-22 


The  final  parameters  of  interest  are  calculated  as 


6f-  <||>  [1  + 8°<*oc  + *c0>] 

(A-52) 

f - ^ 

T 2*ir,T1 

(A-53) 

(A-54) 

After  output  of  IB,  Ic>  h^,  B^,  fT,  VBE,  and  any  other  user-requested 
values,  the  program  control  is  returned  to  the  beginning  of  the  cal- 
culations, and  the  calculations  are  performed  again  with  a new  value  of 


A. 4 Detailed  Subroutine  Structure  - GRTM 

A. 4.1  Module  Usage  and  Interface 

GRTM  is  a general  root  finder  routine  which  is  called  three 
times  during  a pass  through  the  EQSOLVE  equations.  It  is  called  to 
determine  the  onset  parameter  Z,  V',  and  Ic.  Linkage  into  GRTM  is 
accomplished  by  providing  auxiliary  subroutines  which  contain  an 
equation  that  evaluates  to  zero  for  the  proper  value  of  the  variable 
to  be  found.  The  auxiliary  subroutine  names  and  equations  used  in 
EQSOLVE  calls  are  as  follows: 


Name 

Equation 

ZFIND 

(A-2) 

VEFIND 

(A- 5) 

ICFIND 

(A-17) 

A. 4. 2 Subroutine  Input 

Input  to  GRTM  is  the  auxiliary  function  name,  an  initial  guess 
of  roots  requested,  and  an  iteration  convergence  value  described  in 
section  4. 1.4. 4.  The  maximum  iteration  value  is  set  at  50  and  a conver- 
gence criteria  of  1 x 10-^  is  used  for  EQSOLVE. 

A. 4. 3 Subroutine  Output 

The  onset  parameter  Z,  intrinsic  base-emitter  voltage  V_*,  and 

£ 

collector  current  Ic>  are  outputs  of  GRTM  when  called  from  EQSOLVE. 

A. 4. 4 Subroutine  Calculations 

GRTM  uses  Muller's  method  to  find  the  roots  of  f(x).  - This  method 
is  detailed  in  Section  4.3.2. 
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APPENDIX  B 

SENSITIVITY  CURVE  SETS 


B.O 

Introduction 

B.l 

IB  vs- 

VB£  Curves 

B.2 

lQ  vs. 

VBE  Curves 

B.3 

^ vs 

. Ij,  Curves 

B.  A 

Ft  vs. 

1^  Curves 

B.5 

Bf  vs. 

Ic  Curves 

B-l 


B.O  INTRODUCTION 

The  sets  of  parametric  curves  comprising  this  appendix  are 


Calcomp  output  from  the  EQSOLVE  software  code  documented  in  Appendix  A. 

The  curve  sets  show  the  change  in  functional  relationship  of  two-parameters 
of  the  Choma  LSBJT  model  of  a 2 GHz  transistor  under  perturbation  in  value 
of  a third  parameter.  The  five  sets  of  curves  are  selected,  with  the 
exception  of  the  B^  vs.  I set,  as  useful  design-type  curves  which  provide 
the  essential  device  performance  characteristics  at  a glance.  Hence,  the 
effects  of  the  "third"  parameter  variation  on  these  characteristic  curves 
are  produced  by  varying  the  third  parameter  above  and  below  the  nominal 
value.  The  B,  vs.  I curves  are  included  to  provide  insight  into  the 
effects  of  base  pushout  parameter  values  on  the  complex  mathematical 
model  of  base  pushout,  an  effect  which  is  vital  to  understand  for  precise, 
large  signal,  high  injection  modeling. 

The  design  curves  are  included  as  a last  step,  or  "fine-tuning" 
capability  for  a user  who  has  dettermined  a set  of  model  parameters  for  a 
specific  device.  Comparison  of  device  test  data  curves  or  other  CAD 
values  to  mathematical  performance  curves  generated  by  EQSOLVE  or  any 
general  modeling  software  package  (ECAP,  SCEPTRE,  SPICE,  etc.)  will  produce 
curve  differences  similar  to  the  curve  differences  on  one  or  several  of  the 
curve  sets  contained  in  this  appendix.  The  user  then  has  insight  on  the 
erroneous  parameter(s) . 
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ABSTRACT 

This  paper  addresses  the  transient  recovery  characteristics  o£  a 
bipolar  junction  transistor  subjected  to  a saturating  dose-rate  radiation 
event.  An  estimate  of  the  maximum  recovery  time  is  derived  in  terms  of 
parameters  typically  listed  on  device  specification  sheets,  physical  device 
characteristics,  and  circuit  parameters.  Performance  degradation  due  to 
high-injection  phenomena  is  also  addressed. 

1.0  INTRODUCTION 

Saturation  is  invariably  incurred  by  a bipolar  junction  transistor 
exposed  to  an  ionizing  dose  rate  event,  particularly  if  the  transistor 
in  question  conducts  large  quiescent  collector  current  or  if  the  collector 
of  the  device  is  terminated  in  a large  resistance.  Since  a saturated 
transistor  behaves  as  a short-circuit  in  the  sense  that  it  is  incapable 
of  responding  to  input  signal  excitation,  the  time  required  by  the  device 
to  recover  from  the  saturating  excitation  is  clearly  a critical  measure 
of  circuit  performance.  For  example,  this  recovery  certainly  influences 
the  frequency  response  of  a linear  system  and  in  digital  circuitry,  it 
degrades  the  nominal  propagation  delay.  If  the  transistor  in  question  is 
in  a power  supply,  excessive  recovery  requirements  may  well  lead  to 
catastrophic  circuit  or  system  failure. 

This  paper  addresses  the  problem  of  generating  realistic  worst 
case  estimates  of  recovery  time  in  saturated  transistors.  The  equations 
are  cast  as  functions  of  device  and  circuit  parameters.  Although  the 
results  are  predicated  on  simple  modeling  techniques , a semi-quantitative 
account  of  complex  high-injection  phenomena  is  offered  to  refine  the 
first-order  results  obtained. 
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2.0  PHOTOCURRENT  RESPONSE 


Figure  (1)  depicts  a simplified  schematic  diagram  of  a bipolar 

junction  transistor  (BJT)  subjected  to  a gamma  dose  rate  or  gamma  dot 

(y)  event.  The  gamma  dot  event  is  modeled  electrically  by  a collector- 

to-base  photocurrent  source,  i (t).  As  shown  in  Figure  (2),  photocurrent 

PP 

i (t)  is  a simple  current  pulse  of  amplitude  I-,,  and  width  T_ . 

PP  * 

Let  it  be  assumed  that  for  all  steady-state  conditions  corresponding 

to  i (t)  « 0,  the  effective  source  current,  I„,  is  of  a magnitude  appropriate 
PP  * 

to  operation  of  the  BJT  in  the  linear  active  region  of  its  static  charac- 
teristic curves.  Thus,  for  time  t <-Tp,  the  collector  current  is  ic(t)  * ICA* 
The  corresponding  base  current,  1^(0,  is 


ICA/hFE’ 


where  h_  is  the  static  forward  current  transfer  ratio  of  the  bipolar  device. 
r £ 

At  t ■ -Tp,  the  collector  current,  ic(t),  rises  toward  its  saturated 
value,  say  I^g,  in  response  to  abrupt  photocurrent  excitation.  By  inspec- 
tion of  Figure  (1) , this  saturated  value  of  collector  current  is 


V - V 
CC  CSAT 


ipp. 


where  VCSAT  is  the  collector-to-emitter  saturation  voltage.  The  corresponding 
base  current  is  not  linearly  related  to  I^g  by  an  expression  of  the  form 
of  (1)  since  h^g  is  defineable  only  for  the  case  of  transistor  operation 
in  its  linear  active  region.  Accordingly,  if  constant  collector  current 
in  the  amount  of  I^g  is  achieved  in  the  time  frame,  -Tp  < t < 0,  the  base 
current,  IBg,  is 


“ XS  + XPP 


In  (3),  VggAT  is  the  value  of  base-emitter  voltage,  VgE(t),  under  steady- 
state  saturated  conditions. 
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Assuming  chat  Che  transistor  in  Figure  (1)  is  Indeed  saturated 
in  the  neighborhood  of  t < 0,  both  device  junctions  are  forward  biased. 
Resultantly,  a substantial  amount  of  charge  is  stored  in  the  base.  The 
presence  of  excess  base  charge,  which  is  viewed  herewith  as  the  difference 
between  total  base  charge  and  the  charge  that  accrues  if  only  the  base- 
emitter  junction  is  forward  biased,  precludes  an  immediate  response  of  the 
collector  current  to  the  instantaneous  termination  of  the  photocurrent  pulse. 
The  time  required  to  reduce  the  excess  charge  to  zero  or  equivalently, 
the  time  required  to  re-establish  non-forward  bias  across  the  base-collector 
junction,  is  termed  the  "radiation  storage  time,"  tgR. 

The  concept  of  radiation  storage  time  is  best  appreciated  by  an 
investigation  of  Figure  (3),  which  conceptually  plots  the  distribution 
of  charge  injected  into  the  base  under  varying  bias  conditions.  Figure  (3a) 
pertains  to  a BJT  operating  in  its  linear  active  region.  Charge  Qbq,  which 
is  symbolized  by  the  area  of  the  cross-hatched  region,  is  the  equilibrium 
minority  carrier  distribution  evidenced  in  the  base  under  zero  bias 
conditions.  If  the  base  width  (W)  is  small,  the  charge  profile  arising  out 
of  forward  injection  across  the  base-emitter  junction  can  be  approximated 
by  a straight  line.  This  line  has  negative  slope,  since  the  reverse 
bias  existing  across  the  base-collector  junction  promotes  removal  of  charges 
from  base-to-collector.  To  the  extent  that  carrier  transport  resulting 
from  removal  of  base  charge  can  be  explained  in  terms  of  diffusion  mechanics 
alone,  the  collector  current  is  directly  proportional  to  the  magnitude 
of  the  slope  of  Injected  charge  profile.  Thus,  collector  current  in- 
creases are  in  one-to-one  correspondence  with  Increasing  magnitude  of 
charge  profile  slope. 

Under  saturated  conditions,  the  collector  current  is  fixed  at  the 
level  defined  by  (2),  despite  the  fact  that  photocurrent  flow  gives  rise 
to  continuing  injection  of  charge  into  the  base  region.  Since  collector 
current,  and  therefore  charge  profile  slope,  is  constant,  the  resultant 
charge  distribution  must  mirror  Figure  (3b).  The  excess  charge  alluded 
to  earlier  is  Qfig,  which  materializes  primarily  from  injection  across 
a forward  biased  base-collector  junction.  In  short,  both  junctions  are 
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Figure  (3).  Charge  Profiles  In  Base  Region 

(a)  Normal  Active  Operation 

(b)  Saturated  Conditions 

(C)  Threshold  of  Saturated-Normal  Active  Operation 
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forward  biased  In  saturated  regimes,  and  this  condition  is  properly 
reflected  in  Figure  (3b)  by  the  fact  that  the  charge  distribution  is  not 
in  equilibrium  at  either  junction. 

Upon  termination  of  the  radiation  event,  the  collector  current 
remains  nominally  constant  for  a time,  tgR.  This  situation  can  be 
visualized  as  spacially  uniform  removal  of  QgS,  as  suggested  in  Figure  (3c) 
When  tSR  seconds  have  elapsed,  the  collector  current  commences  its  response 
to  photocurrent  pulse  discontinuation,  since  further  removal  of  charge 
mandates  a decrease  in  profile  slope  magnitude.  Note  that  at  the  instant 
when  Qoc  - 0,  the  BJT  is  at  the  threshold  of  re-entering  its  linear 

DO 

active  operating  regime. 

Figure  (4)  offers  plots  of  the  collector  and  base  current  responses 
to  the  photocurrent  excitation  defined  in  Figure  (2).  Observe  that  the 
base  current  decreases  almost  instantaneously  at  t * 0 to  a level,  Igg, 
despite  the  fact  that  the  collector  current  remains  constant  in  the  time 
interval,  0 < t < tgR.  The  magnitude  of  this  change  in  base  current  is 
a function  only  of  Ipp  and  resistive  loading  in  the  base  circuit,  since 
the  voltage  across  the  strongly  charged  base-emittar  junction  cannot 
change  instantaneously.  To  the  extent  that  most  of  the  excess  charge, 

Q , is  swept  primarily  into  either  emitter  or  collector  regions,  the  base 
current  level  of  Ig^  is  essentially  a constant  up  to  time  tgR.  This 
statement  reflects  the  tacit  assumptions  that  (1)  removal  of  Qpg  does 
not  perturb  the  charge  stored  as  a result  of  forward  injection  from 
the  emitter  and  (2)  the  base  current  component  attributed  to  a forward 
biased  base— collector  junction  is  much  smaller  than  the  base  current 
attributed  solely  to  a forward-biased  base-emitter  junction. 

Time  t„  is  termed  the  "recovery  time"  of  the  BJT.  It  is  con- 
ventionally  defined  as  the  time  required  for  collector  current  decrease 
to  within  5%  of  quiescent  value,  following  a saturating  event. 
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3.0  STORAGE  TIME  ANALYSIS 


The  low  Injection  static  equations  which  characterize  BJT 
performance  in  all  operating  regimes  are 


|eVE/VT  -1 

- xs 

V'/V 

e C'  T -1 

hFE 

cVVT-l] 

[eVE/VT  -1 

‘ h 

V' /V 

6 (T  T -1 

_h_' 

-l] 

(3) 

(4) 


where  VC,  and  Vl  are  inttinsic  voltages  across  base-emitter  and  base- 

Ci  L 

collector  junctions,  respectively,  VT  is  the  voltage  equivalent  of 
absolute  junction  temperature  (25.9mV  at  300*K),  and  Ic  is  termed  a reference 
transport  current.  The  static  values  of  emitter  current  IE  and  collector 
current  1^,  flow  in  the  polarity  directions  defined  in  Figure  (1),  and 
voltages  Vl  and  V'  are  positive  when  the  respective  junctions  are  forward 
biased.  Finally,  parameter  h^g  is  the  forward  gain  exploited  in  (1), 
while  h^  is  the  static  reverse  current  transfer  ratio;  i.e.,  h^  is 
the  gain  that  materializes  if  the  roles  of  emitter  and  collector  are 
interchanged. 

An  inspection  of  (3)  and  (4)  readily  indicates  that  both  emitter 
and  collector  currents  are  superpositions  of  the  effects  of  charge  injection 
across  base-emitter  and  base-collector  junctions.  Since  the  terms 
involving  are  dominant  for  forward  or  linear  active  modes  of  operation, 
while  the  terms  in  7'  predominate  under  inverted  operational  circumstancas, 

w 

it  is  reasonable  to  decompose  Ig  and  1^.  into  forward  and  inverse  current 
components.  Accordingly,  let 


XE  " ^F  " ^I 

I m I - I 

AC  aCF  ^CI’ 


(5) 

(6) 
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In  (7)  through  (10), 


XCI  " 'W'V 


hFE 

^Ve  + 1 
bFR 

hFR  + 1 


(10) 

(ID 

(12) 


The  decomposition  of  (3)  and  (4)  into  the  form  of  (5)  and  (6)  is 
significant  in  the  sense  that  the  constituents  of  the  BJT  current  can  be 
identified  in  terms  of  the  charge  storage  components  which  materialize 
in  the  base  due  to  injection  across  either  junction.  In  particular, 
if  Ogg  in  Figure  (3)  can  be  attributed  solely  to  charges  which  traverse 
a forward-biased  (V£  > 0)  base-collector  junction,  both  IRj  and  1^  are 
non-zero  only  if  the  time  rate  of  change  of  QRg  is  non-zero.  It  follows 
that  upon  termination  of  the  photocurrent  pulse,  a solution  for  radiation 
storage  time,  tgR,  embodies  an  investigation  of  the  time  at  which  IRj 
and  Ig^  vanish. 

Although  the  foregoing  solution  tact  is  meaningful,  (5)  through 
(10)  are  Incapable  of  delivering  time-domain  response  information  unless 
they  are  modified  to  account  for  charge  dynamics  in  the  base.  These 
dynamics  are  typically  modeled  by  currents  flowing  across  nonlinear 
junction  capacitances  which  charge  or  discharge  upon  application  or 
removal  of  the  photocurrent  pulse.  A first  order  account  of  dynamical 
charge  effects  is  accomplished  by  replacing  constants  a R and  aR  in  (11) 
and  (12)  by  the  respective  complex  frequency  transforms. 
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(13) 


jr. 


i 


°f<8>  - TT77Z; 


(14) 


I 


In  the  above  relationships.  Is  defined  as  the  forward  galn-bandwldth 

product,  while  u>R  is  the  inverted  analog  of  Specifically,  uF  is 

the  radial  gain  bandwidth  product  of  forward  common-emitter  current 

gain,  and  u>0  is  the  radial  gain  bandwidth  product  of  Inverse  common- 
K 

emitter  current  gain. 

In  the  time  interval,  0 <_  t <_  tgR,  it  may  be  argued  that  the 
inverse  components  of  emitter  and  collector  currents  are  expressible 
in  the  form, 

hlU)  “ + AiEl(t> 

iCI(t)  " ici(0)  + Aici(t)  * 


(15) 

(16) 


In  (15)  and  (16),  iEI(0)  and  iCI(0)  are  steady_8tate  inverse  currents 
which  exist  immediately  prior  to  termination  of  photocurrent  excitation. 
On  the  other  hand,  Ai£I(t)  and  AiCI(t)  are  changes  incurred  in  the 
inverse  emitter  and  collector  currents  during  the  interval  of  radiation 
storage  time.  Clearly,  radiation  storage  time,  tgR,  can  be  defined 
implicitly  as  the  largest  of  the  two  times,  tgpp.  and  tgR(,,  such  that 


AiEI^tSRE^  " *11^' 


(17) 


AiCI(tSRC*  " *CI*0** 


(18) 


It  is  necessary  to  stipulate  \ 

CSR  " MAX^tSRE’  tSRC^’ 

since  the  transistor  re-enters  its  active  operational  region  only  when 
all  Inverse  current  flow  ceases. 
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From  (5)  through  (10) , it  is  easily  Shown  that 


Vs)  f 

hl(a)  m T-  aF(s)Vs)  fF(s)IE 


ap(s)IE(s)  - Ic(s) 


*«<"> 


[S 


(s)IE(s)  - Ic(s) 
- ap(s)aR(s) 


whence,  by  (13)  and  (14) 


wd**  - 


( X T — T 

i (0)  = I - FQ  E lC 

1CIQ0;  " XCI  l - a a * (23) 

FO  RO 

In  (22)  and  (23),  it  is  to  be  understood  that  I£  and  Ic  represent  the 
steady  state  irradiated  values  of  emitter  and  collector  current,  respectively. 

Since  the  collector  current  is  nominally  constant  for  0 < t < tgR, 
the  transforms  of  the  incurred  changes  in  inverse  currents  are,  from 
(20)  and  (21), 

a (s)a  (s) 

SIEI<»>  ■ rr  <24> 


ilcl(s)  - AIEI(.)/aR(s). 

From  Figures  (1)  and  (5),  the  net  change  in  emitter  current  during  the 
storage  time  interval  is 

' * [(*“  + ' (rBS  + Xcs)]-  -7  hs  - fas  - 1Bo)]  • 
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in  the  interval,  0 <_  t <_  tgR(,,  for  the  change  in  the  inverse  component  of 


collector  current. 


(1)  The  inverse  bandwidth  is  much  smaller  than  the  forward  bandwidth 
owing  to  the  reduced  transport  efficiency  incurred  as  a result  of 
charge  injection  from  a lightly  doped  collector  and  charge  collection 
by  a heavily  doped  emitter. 


n 1 n 


An  expression  for  tg^  evolves  straight-forwardly  if  (17),  (22), 
and  (31)  are  combined.  A satisfying  form  for  this  expression  rests  on 
the  observation  that  in  (22), 


“FC^E  XC  “ aFO(1CS  + ~ *CS  ” aFO(IBS  ” XCS ^hFE ^ * 


In  particular,  it  can  be  shown  that 


tSRE 


» 


with 


Similarly, 


A ZBS  ~ TCS/hFE 
0 IBS  “ 1B0 


Since  a-poa-RQ  is  less  than  unity»  tgRE  > tSRC  ^ thus. 


tSR  5 fcSRE* 


(33) 


(34) 


(35) 


(36) 


It  might  be  noted  in  passing  that  the  Inequality,  tgRE  > tgRC  is  a 
reasonable  result.  This  assertion  follows  from  the  fact  that  since 
Inverse  current  flow  reflects  charge  transport  from  collector-to-emitter, 
the  time  difference,  tg^  - tgR(,,  represents  delays  associated  with 
charge  transport  across  the  base  layer  of  the  BJT. 

There  is  considerable  engineering  significance  to  the  parameter, 

N , defined  by  (34).  Observe  that  the  ratio,  ^S^FE*  is  the  minimun 
base  current  required  for  BJT  saturation  in  the  steady  state.  This  is 
to  say  that  a static  base  current  of  Is  commensurate  with  the 

onset  of  Inverse  injection  across  base-collector  junction.  A slightly 
larger  base  drive  delivers  QRg  > 0,  and  the  BJT  is  forced  into  saturation. 
It  follows  that  since  IRg  is  the  actual  steady  state  base  current  for 
the  saturated  BJT,  (IBS  - 1Cs/hyE^  in  the  numerator  on  the  right  hand  side 
of  (34)  represents  the  excess  base  current  commensurate  with  assurance 


of  a saturated  condition.  Recalling  that  I^q  is  the  base  current 
evidenced  Immediately  after  cessation  of  the  photocurrent,  (Igg  - I^q) 
in  (34)  can  be  viewed  as  the  driving  current  serving  to  effect  recovery 
of  the  BJT.  Since  IBQ  must  be  smaller  than  the  BJT  is  indeed 

to  recover,  Nq  is  always  less  than  one.  Accordingly,  Nq  might  be  termed 
a "recovery"  factor  in  the  sense  that  it  compares  the  base  current 
needed  for  device  recovery  to  the  excess  base  current  required  to  ensure 
device  saturation.  As  is  inferred  by  Figure  (5),  minimal  radiation 
storage  time  is  in  one-to-one  correspondence  with  small  Nq.  In  turn, 
small  Nq  demands  a small  Rg/Rgg  ratio,  since  from  (3),  (27),  and  (34), 

/ fjL-Vi  VBSAT/RBB  * ICS/hFE  ~ IS 

°m\  vi 


(37) 
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RECOVERY  TIME 


From  Figure  (4),  the  recovery  time,  tR,  is  the  superposition  of 
radiation  storage  time,  tgR,  and  the  so-called  fall  time,  tR.  In  the 
fall-time  interval,  tgR  < t < tR,  the  BJT  is  in  its  linear  active  region 
of  operation  and  if  pole  dominance  can  be  assumed,  the  collector  -current 
waveform  may  be  approximated  by  a single  time  constant  expression  of  the 
form. 


ic(t> 


*CA  - (ICA  - W''  ‘ 'SR)/Tr 


u(t  " tSR)' 


(38) 


In  (38),  y(t)  is  the  unit  step  function,  and  TR  is  the  first  time  moment 
of  the  small-signal  model  used  to  simulate  the  collector  current  response. 
Taking  ic(tR)  * 1.051^,  (38)  delivers. 


t 


F 


(39) 


Figure  (6)  depicts  the  small-signal  model  alluded  to  in  the 
preceeding  paragraph  [1].  In  this  model,  r is  the  diffusion  resistance 

IT 

of  the  base-emitter  junction,  tq  is  the  incremental  output  resistance,  Cy 
is  the  transition  capacitance  of  the  reverse-biased  base-collector  junction, 
is  the  diffusion  capacity  of  the  base-emitter  junction,  and  SQ  is 
the  small-signal  common- emit ter  current  gain  of  the  BJT.  The  first 
time  moment,  which  equates  to  the  inverse  of  the  dominant  pole  frequency, 
is  the  sum  of  the  time  constants  individually  attributed  to  and  to 
C^,  respectively.  By  straight-forward  circuit  analysis, 

+ cw]  + t40) 


where 


**+RB  + RBB  * 


(41) 


“l 


roRC 
ro  + RC 


(42) 
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Assuming  6q  » r^/R^,  (40)  can  be  cast  into  the  form 


Ip  -O  ****  r»  ) (“F  + “lC“)  6°' 


where 


F - r (c  + C ) 

IT  17  V 

is  the  short-circuit  common-emitter  gain-bandwidth  product,  as  defined 
in  conjunction  with  (13) . 

Care  must  be  exercised  in  the  utilization  of  (43)  and  (44) , since 
these  results  derive  from  a linearized  model  which  may  inappropriately 
simulate  the  effects  of  large  collector  current  excursions.  In  particular, 
the  model  does  not  incorporate  the  variations  in  small-signal  parameters 
incurred  by  large  changes  in  device  currents  and  voltages.  Accordingly, 
one  can  rationalize  the  tact  of  conservatively  estimating  tp  by  choosing 
a worst  case  maximum  value  of  Tp.  From  (43)  and  its  companion  equations, 
it  is  clear  that  a realistic  maximum  value  for  Tp  is 


"““““[•W  +RcC«>']’ 


with  the  understanding  that  BqMAX  and  “fmIN  are  extrema  Pertinent  to 
the  current  excursion,  I^-to-I^,  and  the  corresponding  collector- 
emitter  voltage  variation,  ^CSAT“to“ ^VCC  “ ICARC^*  Furthernore » c0t,» 
the  common-base  output  capacitance  for  zero  collector-base  bias,  is 
a commonly  specified  device  parameter  which  exceeds,  but  closely 
approximates,  small-signal  capacitance  C . 
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EXAMPLE  CALCULATION 


Assume  that  the  circuit  of  Figure  (1)  has  Rc  ■ 500  ohms, 

Rjg  * 5,000  ohms,  Vcc  - 5 volts,  and  Ig  “ 200  micro-amperes . Assume 
further  that  the  BJT  is  characterized  by  the  following  parameters: 


byE 

fTMIN 

fR 

Cob 

0oMAX 


bFR 

V 

BSAT 

V 

CSAT 


50; 

80MHz 

3MHz  (inverse  gain-bandwidth) ; 

2.5pF; 

150; 

800  ohms 

0.4; 

2.3  volts  (includes  ohmic  drop  in  Rg) ; 
0.55  volts. 


Finally,  let  the  BJT  be  exposed  to  a gamma  dose- rate  event  which  gives 
rise  to  a photocurrent  pulse  amplitude  of  2 milll-amperes . 

From  (2) , the  collector  saturation  current  is 


rCS  - 5 500,S5  - «><10'3>  ■ *■»*' 


while  (3)  yields  for  saturated  base  current, 

\ 

IBS  - (200)  (10'6)  + (2)  (10~3)  - - 1.74mA. 

Talcing  V__  - 0.78  volt  for  linear  active  operation.  Figure  (1)  shows 
BE 

that  the  pre-irradiated  value  of  base  current  is 


D-20 


<v 


hi.  * h " i?  ‘ «°0>  <10'6>  - 3555  ‘ 44uAi 


hence,  from  (1),  the  corresponding  collector  current  is 


XCA  " (50>(44>  “ 2.2mA. 


Finally,  (27)  produces 


» (1.74)  (10“3)  - (|iSS)  (2)  (10-' J)  - 15.9yA 


-3, 


BO 


5800 


Using  (11)  and  (12) , 


and 


Then  by  (30) , 


“po  - i - °-96 


“so  ■ ■ °-29- 


- |l  - (0.98)  (0.29>  1 1 2,  (3)  (106) j - 13.5MRPS. 


Equation  (34)  gives  for  the  the  recovery  factor, 


N „ (1. 74) (10~3)  - (6.9) (10~3)/50  _ 9 g 

* (1.74) (10-3)  - (IS. 9) (10‘6)  ’ 


whence  (33)  and  (36)  deliver 


'SR 


(13.5)  do6)  U I1'  °'929l 


- 196nsec 


for  the  estimated  radiation  storage  time. 


Using  (21), 


■ 486nsec. 

Then  from  (39), 

CFMAX  “ (^86) (10~9)  | 3.0  + Ln  - lj  - 1.83vsec. 

Thus,  the  estimated  worst  case  BJT  recovery  time  is 

Wx  ■ ‘nux  + CSE  • 2-°“sec- 

Observe  in  the  foregoing  example  that  the  storage  time  is  in 
the  order  of  10%  of  the  maximum  recovery  time.  This  result  generally 
reflects  situations  actually  encountered  in  realistic  radiation  test 
environments.  Accordingly,  it  follows  that  BJT  recovery  time  is  largely 
determined  by  time  moment  parameter  ^ turn,  is  minimal 

for  devices  which  exude  low  collector-base  capacitance,  large  gain 
bnadwidth  product,  and  low  small-signal  current  gain. 


TFMAX  ” 150 


2ir  (80)  (10°) 


+ (500)  (2.5) (10"12) 


6.0  CONCLUSIONS 


An  analysis  of  the  radiation  recovery  transient  in  a bipolar 
junction  transistor  subjected  to  a saturating  dose  rate  event  has  been 
developed.  The  analysis,  which  is  supplemented  by  exemplifying  cal- 
culations of  recovery  time,  reveals  device  parameters  and  performance 
attributes  that  must  be  critically  scrutinized  if  recovery  time  is  to 
be  constrained  to  a minimum.  Although  the  general  nature  of  sensitive 
device  parameters  is  properly  uncovered,  it  is  crucial  that  approximations 
invoked  during  the  course  of  analysis  be  understood,  particularly  if  the 
primary  objective  of  utilizing  the  analysis  in  a given  application  is 
realistic  absolute  estimation  of  recovery  time. 

Undoubtedly,  the  most  significant  assumption  is  that  BJT  currents 
flow  by  diffusion  mechanisms  alone.  In  simplified  terminology,  it  has 
been  presumed  that  (1)  the  injected  charge  profile  in  the  base  is 
linearly  graded,  (2)  the  collector  current  is  proportional  to  the 
magnitude  of  profile  gradient,  and  (3)  the  excess  stored  charge  is 
removed  prior  to  removal  of  charge  accumulated  from  forward  Injection 
across  the  base-emitter  junction.  In  reality,  drift  components  of  BJT 
currents  are  at  least  comparable  to,  if  not  larger  than,  the  diffusion 
components  in  saturated  devices.  These  drift  components  loom  increasingly 
more  important  if  the  actual  base  current  is  allowed  to  increase  sig- 
nificantly over  the  minimum  base  current  (I^g/h^g)  commensurate  with 
the  threshold  of  saturation. 

Aside  from  compromising  the  foregoing  three  presumptions,  the 
immediate  effect  of  non-negligible  current  transport  via  high-field 
drift  mechanisms  is  high-current,  low-voltage  degradation  in  the  gain- 
bandwidth  product,  . The  maximum  amount  of  degradation  is  propor- 

f 

tlonal  to  the  square  of  the  ratio  of  epitaxial  layer  length- to-neutral 
base  width.  Depending  on  the  type  of  BJT  addressed,  the  extent  of  this 
degradation  can  be  in  the  range  of  lS-to-20.  Similar  attenuations  are 
evidenced  with  respect  to  both  the  static  and  small-signal  current  gains, 
although  the  factor  by  which  these  gains  degrade  is  not  as  severe  as  the 
degradations  quoted  for  gain-bandwidth  product.  The  upshot  of  the  matter 
is  that  if  manufacturer's  gain-bandwidth  product  and  gain  specifications 


r 


are  used  to  compute  T_,.,,,  it  is  conceivable  that  the  estimated  maximum 
recovery  time  might  be  low  by  factors  in  the  range  of  5-to-10,  since 
device  specifications  are  invariably  applicable  to  operational  regions 
other  than  high- injection  regimes. 

It  must  also  be  remembered  that  (38)  and  (39)  presume  device 
linearity  and  pole  dominance  for  tgR  < t < tR.  The  effects  of  inherent 
device  nonlinearities  are  especially  pronounced  if  the  saturated-to- 
quiescent  collector  current  ratio,  is  larger  than  four  or  five. 

Ringing  in  the  collector  and  base  current  responses  is  indicative  of 
either  circuit  nonlinearity  or  a non-dominant  pole  population.  Indeed, 
if  non-linearity  is  particularly  severe,  the  analytical  concept  of  a 
pole  cannot  be  exploited  in  the  computation  of  fall  time  tR.  It  can  be 
envisaged  that  the  impropriety  of  either  linearity  or  pole  dominance 
is  conducive  to  recovery  time  estimates  that  are  low  by  20Z-to-40%. 

In  summary,  it  is  suggested  that  the  foregoing  equations  are 

useful  in  ascertaining  device  characteristics  which  strongly  influence 

» 

the  radiation  recovery  response  of  a BJT.  They  are  also  useful  when 
applied  to  the  problem  of  bracketing  a realistic  estimate  of  maximum 
recovery  time.  The  lowest  anticipated  maximum  recovery  time  derives  from 
direct  employment  of  manufacturer's  nominal  specifications  for  key  device 
parameters.  The  largest  estimate  for  maximal  recovery  time  is  obtained 
by  allowing  gain  and  gain-bandwidth  product  to  degrade  in  accordance 
with  specified,  or  estimated  high-inj ection  device  performance 
characteristics. 
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ABSTRACT: 


This  paper  proposes  the  use  of  an  existing  large-signal  circuits 
model  to  simulate  bipolar  transistor  performance  degradation  induced  by 
neutron  irradiation.  It  is  shown  that  only  three  radiation  parameters  are 
required  to  establish  a model  for  neutron  effects  over  the  entire  practical 
current  and  voltage  ranges  of  bipolar  transistor  operation. 
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The  predominant  effects  of  neutron  displacement  damage  In  silicon  are 
an  Increase  In  carrier  generation  rate  and  a decrease  In  majority  carrier 
concentration.  In  a bipolar  junction  transistor  (BJT) , these  effects  combine 
to  produce  current  gain  attenuation,  degradation  In  gain-bandwidth  product, 
enhanced  base-collector  leakage  current,  increased  collector  breakdown  voltage, 
and  Increased  collector  saturation  voltage  [1].  Additionally,  neutron  dis- 
placement damage  in  a BJT  may  significantly  perturb  the  nominal  thermal 
characteristics  of  a bipolar  circuit.  Depending  on  circuit  configuration, 
the  nominal  thermal  behavior  can  either  be  aggravated  or  Improved  by  neutron 
bombardment  [2]- [3]. 

Traditionally,  the  performance  of  bipolar  circuits  immersed  in  a 
neutron  environment  is  simulated  by  exploiting  various  damage  constants  which 
are  enumerated  from  experimental  device  characterizations.  In  particular, 
gain,  leakage  current,  and  gain-bandwidth  product  are  measured  prior  to  and 
immediately  following  controlled  neutron  irradiation  of  a given  bipolar  device. 
This  data  is  used  to  extrapolate  simple  relationships  between  pre-irradiated 
and  irradiated  values  of  various  device  parameters.  An  example  of  such  a 
relationship  is  the  Messenger-Spratt  equation  [4], 


£ +Vfo*> 


(1) 


where  8.  is  the  irradiated  value  of  static  short-circuit  current  transfer  ratio, 
9 

8 . t is  the  transit  time  of  minority  carriers  across  the  neutral  base,  9 is 

0 r U 

neutron  fluence,  and  is  the  pertinent  damage  constant.  Other  device  per- 
formance indices,  such  as  gain-bandwidth  product,  breakdown  voltage,  and 
leakage  current,  are  related  to  device  neutron  fluence  exposure  by  expressions 
whose  forms  are  analogous  to  (1) . 

Although  neutron-induced  damage  assessment  equations  are  useful  both  for 
simulating  the  dependence  of  Irradiated  electrical  parameters  on  their  pre- 
irradiated  counterparts  and  for  establishing  a basis  of  neutron  hardness 
screening,  at  least  two  fundamental  objections  can  be  raised  if  these  equations 
are  used  exclusively  in  evaluating  bipolar  response  to  neutron  irradiation. 
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First,  the  damage  constant  methodology  can  only  simulate  measured  results. 

It  cannot  be  applied  to  the  problem  of  predicting  the  radiation  response  of 
bipolar  transistors,  because  the  damage  constants  are  empirical  performance 
barometers.  These  constants  are  not  well-defined  in  terms  of  intrinsic  physical 
or  geometrical  device  properties  but  instead,  they  derive  purely  from  experi- 
mental data. 


A second  shortcoming  is  that  most  of  the  damage  constants  invoked  in  a 
neutron  radiation  response  exercise  are  sensitive  functions  of  device  injection 
level.  For  example,  K.  in  (1)  varies  with  the  voltage  and  current  at  which 
gain  is  measured.  The  upshot  of  the  matter  is  that  most  damage  assessment 

equations  are  applicable  only  to  the  injection  regime  in  which  they  are  for- 
mulated. Attempts  made  to  invoke  these  same  equations  to  damage  assessment 
at  significantly  lower  or  higher  injection  levels  generally  result  in  sub- 
stantial analytical  errors.  It  follows  that  the  methodology  in  question 
conveys  little,  if  any,  understanding  of  the  sensitivity  of  neutron  radiation 
response  with  respect  to  the  possibly  wide  Injection  range  which  a bipolar 
device  encounters  in  a prescribed  circuit  application. 

This  paper  presents  a new  bipolar  neutron  damage  model  which  is 
uniformly  applicable  to  the  low- through-high  injection  range  of  BJT  operation. 
The  model  is  an  extension  of  an  existing  circuits  model  which  is  capable 
of  predicting  large-signal  transient  responses  to  forcing  inputs  that  can 
drive  the  transistor  from  low-to-high  injection  regimes.  In  addition  to  the 
model  parameters  pertinent  to  simulation  and  prediction  of  traditional 
electrical  responses,  only  three  additional  parameters  are  required  for 
an  evaluation  of  the  effects  of  neutron  radiation.  These  parameters  are 
independent  of  voltage  and  current  and  moreover,  they  can  be  related  to 
damage  constants  which  have  been  enumerated  from  a damage  assessment  appro- 
priate to  mid-injection  operational  regimes. 


2.0 


FUNDAMENTAL  TRANSISTOR  MODEL 


The  bipolar  transistor  circuit  model  diagrammed  in  Figure  (1)  is 
proposed  as  the  foundation  for  neutron  radiation  response  studies.  ^ The 
model  possesses  all  of  the  attributes  of  the  classical  Ebers-Moll  [5]  and 
Gummel-Poon  [6]  topologies  and  in  addition,  it  provides  the  following  sim- 
ulation capabilities. 

(1)  . The  effects  of  base-collector  depletion  layer  width  modulations  arising 

out  of  collector  junction  voltage  changes  are  Incorporated . This 
phenomenon,  known  as  the  Early  Effect  [7],  gives  rise  to  current  gain 
dependence  on  collector  voltage.  Analogously,  an  account  is  made  of 
Inverse  Early  Effect,  which  produces  inverse  beta  dependence  on  base- 
emitter  voltage. 

(2) .  Low-injection  forward  current  gain  attenuation  due  to  non-ideal  diode 

currents  flowing  across  the  base-emitter  junction  is  replicated  by  the 
model.  In  Figure  (1),  the  lumped  equivalent  to  the  non-ideal  base- 
emitter  diode  current  is 

IBEB  ■ ier<e'Wie  -U  • <2> 

where  is  junction  thermal  voltage,  and  and  N£  relate  implicitly 
to  emitter  efficiency,  degree  of  emitter  sidewall  injection,  and  the 
recombination  velocity  at  the  base-emitter  interface.  Parameter  N£  is 
always  greater  than  or  equal  to  unity.  The  diode  current,  IBC^,  mathe- 
matically parallels  the  form  of  (2)  and  allows  for  low-level  inverse 
current  gain  attenuation  for  transistors  operated  in  inverted  or 
saturated  modes. 


(1)  The  model,  which  was  developed  by  the  author  at  the  University  of 

Pennsylvania,  is  embedded  in  the  ISPICE  computer-aided  design  program 
and  is  available  on  a timeshare  network  marketed  by  National  CSS,  Inc., 
Norwalk,  Connecticut. 
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(3).  High- injection  forward  current  gain  degradation  caused  by  conductivity 
modulation  in  both  the  base  and  quasi-saturated  collector  regions  is 
incorporated  in  the  model  of  Figure  (1) . The  element  germane  to 
this  simulation  capability  is 


V'/V 

Is(eV  TE  -1), 


(3) 


where  Q_,  the  net  majority  charge  in  the  effective  base,  has  a value 
of  QB0  under  zero  bias  (V£  * ■ 0)  conditions.  The  charge  function, 

Qg,  which  is  the  superposition  of  charge  Q£  stored  in  the  base-emitter 
transition  layer,  charge  stored  in  the  base-collector  depletion 
region,  and  the  charges  (Q^,  and  Q^)  required  to  preserve  space-charge 
neutrality  in  the  face  of  charge  injection  across  either  junction, 
is  given  by 


t I (eVVTE  -1)  x I (eVC/VTC  -1) 

<1  + V+JLL-75 


B 


In  (4), 


(4) 


(5) 


where  V_  and  V are  inverse  and  forward  Early  effect  parameters, 

ER  CR 

and  t (t  ) is  the  forward  (inverse)  transit  time  of  minority  carriers 
E R 

in  the  field-neutral  base.  Because  of  base  pushout  phenomena,  these 
transit  times  vary  with  device  current  and  voltage  [8].  Parameter 
Ig  in  (3)  and  (4)  is  directly  proportional  to  the  product  of  intrinsic 
carrier  concentration  and  carrier  diffusivlty  in  the  base,  and  it 
relates  Inversely  to  the  total  majority  charge  in  the  base  region. 


(4).  High-injection  gain-bandwidth  product  (fT)  degradation  caused  by  base 
pushout  and  conductivity  modulation  in  the  collector  layer  is  a model 
feature.  The  model  also  addresses  the  related  nonlinearity  in  effective 
collector  resistance,  (RC(,  + R^q)  . 
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(5) .  The  effects  of  static  emitter  crowding  on  the  effective  base  resistance, 

(RgQ  + Rgjj)»  are  embodied  in  the  model.  Excess  phase  and  both  static 
and  dynamic  conduction  across  the  collector  substrate  interface  are 
also  included. 

(6) .  A first  order  approximation  of  the  distributed  effects  of  high-frequency 

excitation  is  introduced  by  way  of  the  semi-empirical  parameter, 

3.0  LOW- INJECTION  CURRENT  GAIN  DEGRADATION 

Under  static  forward  bias  operating  conditions,  all  charge  functions 
vanish,  the  substrate  diode  current,  Igg,  is  very  nearly  zero,  and  currents 
Ibcn>  and  IIT*  w^ic^  are  exponential  functions  of  the  intrinsic  base- 

collector  junction  voltage,  are  essentially  zero.  Accordingly,  Figure  (1) 
verifies  that  base  current  Ifi  is  given  by 


*B  " *BE  + IBEN* 

where  1^^  is  defined  by  (2)  and 

Tbe  " J5,  teVE/VTE  -1] . (7) 

0 

In  (7),  Bq  is  the  maximum  value  of  the  static  transfer  ratio, 

- § • <8) 

under  the  special  condition  of  zero  bias  (V*  ■ 0)  across  base-collector  junction. 

Returning  to  Figure  (1) , it  is  clear  that  under  the  conditions 
stipulated  above,  the  collector  current,  Iq,  is  identical  to  the  current, 

1^.  From  (3)  and  (A)  it  can  be  shown  that 
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Is(eVE/VTE  -1) 


1 +X, 


-1  *-X")  + Q^CeV^  -1)  + Qgjj (eVC/VTC  -1) 


1/2  ' 


(9) 


where 
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The  term  involving  QOVI  in  (9)  is  essentially  zero,  provided  that  the  BJT 

KXi 

is  operated  in  its  linear  active  domain.  Moreover,  for  low-to-moderate 
injection  levels,  the  term  containing  is  negligible,  so  that  (9)  collapses 
to 

xc  e^/yiE-  (12> 

In  (12),  it  is  assumed  that  the  exponential  term  is  much  larger  than  unity. 

If  this  same  reasonable  assumption  is  invoked  on  (2) , (8)  can  be  shown  to 
deliver 


with 


(13) 


(14) 
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Figure  (2)  portrays  the  dependence  of  hgg  and  V',  and  hence  the 
dependence  of  hgg  on  the  natural  logarithm  of  collector  current  1^,.  Re- 
calling (5),  the  effect  of  parameter  X£,  as  depicted  in  the  figure,  is  to 
simulate  the  sensitivity  of  current  gain  with  respect  to  the  reverse  biasing 
voltage  applied  at  the  collector.  Observe  then  that  for  increasing  V', 
hpE  correctly  approaches  a constant  at  fixed  collector  junction  voltage, 
provided  that  the  contribution  of  V*  to  XE  is  insignificant  in  comparison 
to  the  contribution  of  V'.  For  decreasing  forward  bias,  hgg  attenuates 
in  such  a way  that  at  V'  • V^,,  hgg  is  one-half  of  the  current  gain  at 
moderate  current  levels.  Notice  that  Vg  is  dependent  on  both  Igg  and  Ng, 
which  implicitly  relate  to  the  characteristics  of  the  base-emitter  interface. 
Finally,  Ng  fixes  the  slope  of  the  gain  curve  at  Vg  ■ Vg,  If  Ng  ■ 1,  which 
infers  only  idealized  diode  current  flow  across  the  base-emitter  junction, 
hgg  is  constant  over  the  addressed  range  of  Vg, 

The  preceding  discussion  applies  to  a BJT  operated  in  a traditional, 
non-nuclear,  electrical  environment.  When  this  environment  is  perturbed 
by  significant  neutron  fluence,  one  can  postulate  that  additional  base  current, 
say  dig,  must  be  supplied  in  order  to  maintain  the  collector  current  at  its 
pre- irradiated  level.  Let  this  base  current  enhancement  be  expressed  as 


AIB  “ % 


|eVE/VTE  -lj 


V'/N  V 

+ I.  e V E TE 

$o 
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where  it  is  understood  that  I.  is  a function  of  the  neutron  energy  imparted 

9 

to  the  active  region  of  the  base-emitter  junction.  On  the  other  hand,  I^Q 
is  parametrically  related  to  the  neutron  energy  absorbed  by  the  emitter 
sidewalls,  the  emitter-base  depletion  region,  and  the  surface  of  the  depleted 
layer  associated  with  the  base-emitter  junction. 

Clearly,  the  induced  increment  in  base  current  can  be  modeled  by 
incorporating  two  additional  current  sources  in  the  base-emitter  circuit 
of  Figure  (1).  Current  IgEN  is  placed  in  parallel  with  a source  , whose 
current  is  defined  by  the  second  term  on  the  right  hand  side  of  (15) , while 
the  first  term  in  this  relationship  is  a diode  in  shunt  with  current  IgE 
in  Figure  (1).  Notice,  however,  that  these  modifications  infer  a mere 
redefinition  of  parameters  Pq  and  Igg.  Specifically,  SQ  and  TER  in  the 
traditional  model  can  be  replaced  by  irradiated  parameter  values,  B^  and 
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$R,  respectively,  such  that 


♦ 1 + Wh 


" rER  + V 


A number  of  noteworthy  comments  can  be  offered  at  this  juncture.  First, 

only  two  constant  parameters,  1.  and  1.  , are  required  to  establish  neutron- 

$ $0 

induced  changes  In  the  current  gain  response  over  the  range  of  very  low-to- 

moderate  Injection  levels.  In  particular,  parameter  1^  serves  to  attenuate 

the  maximum  current  transfer  ratio  at  moderate  current  levels,  while  I. 

. $o 

Influences  the  nature  of  h^  in  low  injection  regimes.  The  latter  point 
becomes  clear  when  one  replaces  1^  and  8q  in  (14)  by  I R and  B to  obtain 


V i 1 * ho/LEBL 

TE  1 + 8o  V^S 


K \ne  - 1 


Thus,  the  base-emitter  junction  bias  commensurate  with  an  hpg  that  is  one- 
half  of  the  maximum  current  transfer  ratio  at  moderate  current  levels  shifts 


as  a function  of  the  radiation  parameters,  1^  and  I^q.  Note,  however,  that 
while  the  maximum  current  transfer  ratio  always  decreases  In  response  to 
neutron  Irradiation,  the  half-gain  junction  bias  level  does  not  necessarily 
change,  since  if  I^/I^  ■ B^/Ig,  = VR.  In  effect,  the  voltage  change, 
- VK,  monitors  the  relative  deposition  of  energy  between  the  active 
and  inactive  portions  of  the  base-emitter  junction. 

A second  interesting  point  is  that  (16)  can  be  written  as 


±-mk-  + 3. 

*9  So  *S  ’ 


which  is  precisely  of  the  form  of  the  classical  Messenger-Spratt  relationship. 

Indeed,  a comparison  of  (1)  and  (19)  infers  the  I,  relates  to  neutron  damage 

9 

constant  K,  and  fluence  * in  accordance  with 
9 


T 


(20) 
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?» 

k: 
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I.  " K.T--I-#. 
<f>  <(>  F0  S 


It  is  critical  to  note  that  since  1^  is  a constant  model  parameter,  r 


FO 


is  necessarily  a constant.  This  is  to  say  that  Is  the  minority  carrier 


transit  time  across  a field-neutral  base.  Since  a field-neutral  base  is 
ensured  at  moderate  currents  only  under  the  condition  of  a strongly  back- 
biased  collector-base  junction,  and  since  is  inversely  proportional 
to  common-emitter  gain  bandwidth  product  fj,  (20)  is  meaningfully  expressed 
as 


K 


\ “ 


2irf( 


«L 


TMAX  , 


(21) 


where  Is  understood  to  be  the  value  of  f^  appropriate  to  linear  active 

operation  of  a BJT  subjected  to  strong  collector-base  back  bias. 


4.0  HIGH- INJECTION  CHARACTERISTICS 

When  the  base-emitter  junction  is  strongly  forward  biased,  the  term 
involving  in  (9)  cannot  be  ignored.  Moreover,  since  Ng  is  larger  than 
unity,  Igg  is  much  greater  than  and  thus,  the  net  base  current  is 

closely  approximated  by  the  right  hand  side  of  (7) . It  follows  that  for 
V'  « 0,  (7)  and  (9)  combine  to  produce 


^FE 

fL 


„ . . x . qfnic 
(1  + XE)  + 


(22) 


Equation  (22)  is  plotted  in  Figure  (3),  wherein  it  is  understood 
that  zero  collector  current  (Ig  ■ 0)  is  in  one-to-one  correspondence 
with  the  moderate  current  levels  associated  with  regions  on  the  far  right 
hand  side  in  the  graph  of  Figure  (2).  Observe  in  (22)  that  the  high- 

injection  current  transfer  ratio  attenuates  to  one-half  of  its  moderate 

\ 

current  value  at  a collector  current,  say  1^,  given  by 

+ xe>Qb0 


IK- 


(23) 
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Figure  (3).  Forward  Current  Transfer  Ratio  At  High- Injection  Levels. 
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In  arriving  at  (23),  (10)  is  utilized  to  express  1^  in  terns  of  base  region 

transit  tine,  Tp.  It  should  be  pointed  out  that  the  actual  dependence  of 

hpg  on  collector  current  is  more  complicated  than  suggested  by  (22)  and 

(23),  since  base  pushout  at  high-injection  levels  Incurs  a dependence  of 

T_  (and  hence,  Q_.)  on  current  I„  and  the  voltage  parameter,  X . 
t FN  U Ft 

The  development  of  a model  to  simulate  neutron- induced  changes  in 

the  high-injection  current  transfer  characteristics  derives  from  two 

observations.  First,  parameter  0Q  attenuates  to  B^,  as  defined  by  (16). 

Thus,  Bq  in  (22)  is  replaced  by  B^,  thereby  incurring  a uniform  downward 

shift  in  the  hpE  curves  portrayed  in  Figure  (3).  Second,  the.majority 

carrier  concentration  in  the  base  decreases  from  to  in  accordance 

with  the  semi- empirical  relationship  [9] 

^BO 

QB$  " 1 + Kg#  » (24) 

where  Kg  can  be  thought  of  as  charge  damage  constant  which  reflects  neutron- 
induced  decreases  in  majority  carrier  concentration.  Upon  replacement  of 
Q„r  in  (10)  by  Q_. , one  concludes  that  the  irradiated  value  of  Q is 


<W*'W1  + V>-  <25> 

which  directly  supplants  pre-irradiated  parameter  in  (22).  Note,  however, 
that  such  a replacement  of  variables  is  equivalent  to  allowing  the  transit 
time  to  assume  an  irradiated  value,  Tp^,  of 


( 


• V1  + V> 


There  are  at  least  two  important  ramifications  of  (26).  First  and 
foremost,  the  equation  allows  for  the  simulation  of  neutron- induced  per- 
turbations in  the  gain-bandwidth  product,  as  well  as  neutron-induced  changes 
in  the  forward  current  transfer  ratio.  This  statement  exploits  the  simple 
fact  that  f<»r  collector  currents  in  excess  of  a critical  value,  say  Iq(V'), 


fT  “ 2irtp  » :C  ” V 
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In  (27),  Tj,  Is  a function  of  both  collector  current  Ig  and  collector  junction 
bias  voltage,  V' . This  functional  relationship,  which  is  portrayed  graphically 
by  Figure  (4),  is  mathematically  defined  by  a so-called  "base-pushout 
function,"  Bp(Ig,  v') * such  that 


TF  " TFOBF^IC*  V 


bf(ic»  va> 

2irW  * 


(28) 


The  current  parameter,  Iq(V^)  which  defines  the  onset  of  BJT  entry  into 
high- injection  operating  regimes,  and  the  pushout  function,  B^(I^,  V^), 
can  be  related  to  such  device  processing  chaxacterlstics  as  epitaxial 
layer  doping,  base  and  epitaxial  layer  geometries,  minority  carrier 
dlffuslvitles,  and  the  like  [8],  [10].  Alternatively,  the  behavior  of 
Tp  in  high- Injection  domains  can  be  modeled  by  empirical  equations  whose 
parameters  derive  conveniently  from  measured  f ^characteristics  [11].  The 
upshot  of  the  matter  is  that  (26)  allows  for  the  simulation  of  f,p  over  both 
moderate  and  high  injection  regimes  since  if  Tp  in  (27)  is  replaced  by 
Tp^,  (26)  through  (28)  infer  an  Irradiated  value  of  f^  that  is  given  by 


2,,TF0a  + V>W  vc>  1 + V ' 


(29) 


£or  *C  * Io(vc>- 

It  is  useful  to  Interject  that  for  Ig  < IQ (V') , (26)  remains  applicable 
to  the  problem  of  simulating  gain-bandwidth  product  characteristics.  For 
collector  currents  not  exceeding  Iq (Vg) , the  base  pushout  function  is  unity 
and 
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ac 

T 


_1 

+ VTECT^IC^ 


(30) 


where  C^,  the  sum  of  the  transition  capacitances  of  both  junctions,  is 
nominally  insensitive  to  neutron  fluence.  Clearly,  the  neutron  irradiated 
value  of  low-current  f^  is  obtained  by  replacing  tpQ  in  (30)  by  Tp^,  as  defined 
by  (26).  Unlike  the  uniform  neutron-incurred  degradation  of  f^  at  high 
collector  currents,  the  effect  of  neutrons  on  low-current  f^  is  seen  to  be 
non-uniform,  owing  to  the  relatively  small  Influence  which  neutron  fluence 
exerts  on  transition  layer  storage  dynamics. 
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Gain-Bandwidth  Product, 
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Figure  (4).  Common-Emitter  Gain- Bandwidth  Product 
as  a Function  of  Voltage  and  Current. 
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A second  ramification  of  (26)  derives  directly  from  the  observation 
that  this  equation  is  predicated  on  the  fundamental  assumption  that  the 
irradiated  value  of  base  region  majority  charge  is  linearly  related  to  its 
pre-irradiated  counterpart  by  an  expression  of  the  form  of  (24) * Since  both 
in  (11)  and  in  (10)  are  inversely  proportional  to  QM,  it  is  clear 
that  the  irradiated  Inverse  transit  time  behaves  in  a manner  that  is  identical 
to  the  forward  transit  time  behavior  prescribed  by  (26). 

5.0  CONCLUSIONS 


The  fundamental  premise  of  this  paper  is  that  an  existing  and  reasonably 
accessible  electrical  model  of  a bipolar  junction  transistor  is  easily  modified 
to  incorporate  the  cognate  effects  of  BJT  exposure  to  a given  neutron  fluence, 
9.  The  foundation  for  the  required  modifications  is  the  fact  that  (1) 
additional  base  current  must  be  supplied  if  the  collector  current  of  an  ir- 
radiated BJT  is  to  be  maintained  at  its  pre-lrradlated  value,  and  (2)  the 
majority  charge  in  the  active  base  decreases  in  proportion  to  the  neutron 
fluence.  The  additional  base  current  and  charge  perturbation  are  stipulated 
by  (15)  and  (24),  respectively,  which  in  turn  define  three  neutron-related 
modeling  parameters;  namely,  1^,  I^q,  and  Kg. 

Parameter  I.  can  be  determined  in  either  of  two  ways.  One  way 

exploits  (19),  which  entails  a comparison  of  the  pre-lrradlated  and  post- 

irradiated  values  of  the  gain  parameter  BQ,  as  defined  in  Figure  (2).  A 

second  technique  is  to  relate  1^  to  a known  current  gain  damage  constant  in 

accordance  with  (20)  or  (21).  Parameter  I.  requires  an  investigation  of 

90 

current  transfer  ratio  evidenced  at  low  current  levels.  In  particular, 
the  base-emitter  voltage  (which  closely  approximates  the  base-emitter 
junction  bias  at  low  currents)  is  measured  prior  and  subsequent  to  neutron 
irradiation  to  deduce  I^q  by  way  of  (18) . Finally, , Kg  is  most  easily 
estimated  by  monitoring  pre-irradiated  and  irradiated  f^  at  fixed  current 
and  voltage  appropriate  to  operation  in  high-injection  regimes.  The  equation 
of  interest  is  (29). 


It  is  important  to  note  that  while  parameters  I^»  I^Q,  and  Kg  are 
constants,  their  correct  utilization  enables  the  simulation  of  the  damaging 
effects  of  neutron  bombardment  over  all  forward  active  regimes  of  BJT 
operation.  Indeed,  parameter  Kg  is  ostensibly  applicable  to  gain-bandwidth 


product  studies  of  bipolar  devices  operated  In  inverted  or  saturated  regimes. 
Although  inverse  BJT  operation  is  not  specifically  addressed  in  this  work, 
parameters  analogous  to  1^  and  can  presumably  be  found  to  simulate 
inverted  current  transfer  ratio  characteristics.  The  only  problem  inherent 
in  this  presumption  is  the  difficulty  encountered  in  monitoring  pre-irradiated 
inverted  performance  without  incurring  permanent  damage  at  the  base-emitter 
interface . 

It  is  believed  that  the  theory  of  neutron  effects  modeling  developed 
in  this  paper  is  applicable  to  a broad  range  of  BJT  simulation  and  prediction 
problems.  For  example , over  a wide  range  of  practical  bipolar  device  operation, 
the  dependence  of  collector  current  on  temperature  at  fixed  bias  voltage  is 
given  by  [12] 

Ic  - ATMe“VGO/VTE,  # (31) 

where  is  the  bandgap  voltage  and  constants  A and  M are  determined  by 
device  fabrication  processes.  Clearly,  M is  critical  to  the  determination 
of  the  thermal  sensitivity  of  collector  current.  It  turns  out  that  published 
data  related  to  the  dependence  of  minority  carrier  mobility  on  base  region 
impurity  concentration  infers  [2]  that  M behaves  approximately  as 


M * LbM,,,  1^)  - V 


where 


and  Mq  are  constants.  The  irradiated  value  of  M evolves  from 


replacement  of  by  Q^,  whence  by  (31),  the  neutron-induced  thermal 
sensitivity  of  collector  current  is  easily  discerned. 


E-19 


6.0  REFERENCES 

1. 


F.  Larin,  Radiation  Effects  In  Semiconductor  Devices,  New  York: 
John  Wiley  & Sons,  Inc.,  1968,  Chapter  9. 


2. 


S.  M.  Sze,  Physics  of  Semiconductor  Devices,  New  York:  John  Wiley  & 
Sons,  Inc. 


3. 


J.  S.  Brugler,  "Silicon  Transistor  Biasing  for  Linear  Collector 
Current  Temperature  Dependence,"  IEEE  Journal  of  Solid-State 
Circuits,  Vol.  SC-2,  pp.  57-58,  June  1967. 


4. 


G.  Messenger  and  J.  Spratt,  "The  Effects  of  Neutron  Irradiation 
on  Germanium  and  Silicon,"  Proceedings  of  IRE,  pp.  1038-1044, 
June  1958. 


5. 


J.  J.  Ebers  and  J.  L.  Moll,  "Large-Signal  Behavior  of  Junction 
Transistors,"  Proceedings  of  IRE,  Vol.  42,  pp.  1761-1772, 
December  1954. 


6. 


H.  K.  Gummel  and  H.  C.  Poon,  "An  Integral  Charge-Control  Model  of 
Bipolar  Transistors,"  Bell  System  Technical  Journal,  Vol.  49, 
pp.  115-120,  May- June  1970. 


7. 


J.  M.  Early,  "Effects  of  Space-Charge  Layer  Widening  In  Junction 
Transistors,"  Proceedings  of  IRE,  Vol.  40,  pp.  1401-1406, 
November  1952. 


8. 


J.  Choma,  Jr.,  "A  Process-Oriented  Model  for  the  Simulation  of 
Base  Pushout  In  Integrated  Bipolar  Devices,"  IEEE  Transactions 
on  Electron  Devices,  Vol.  ED- 22,  pp.  1079-1086,  December  1975. 


9. 

10. 


Reference  #1,  pp.  157-169. 


J.  Choma,  Jr.,  "Process-Oriented,  High- Injection  Circuit  Models 
for  Integrated  Bipolar  Junction  Transistors,"  ONR  Final  Report, 
No.  N00014-75-C-1177,  November  1976. 


11. 


J.  Choma,  Jr. , "A  Curve-Fitted  Circuits  Model  for  Bipolar  Transistor 
fT  Roll-Off  at  High  Injection  Levels,"  IEEE  Journal  of  Solid-State 
Circuits,  Vol.  SC-11,  pp.  346-348,  April  1976. 


12. 


K.  E.  Kuijk,  "A  Precision  Reference  Voltage  Source,"  IEEE  Journal 
of  Solid-State  Circuits.  Vol.  SC-8,  pp.  222-225,  June  1973. 


E-20 


APPENDIX  F 

ACKNOWLEDGEMENTS 


F-l 


= 


The  number  of  people  who  contributed  to  this  research  project 
are  too  numerous  to  list  herewith.  However,  I would  be  neglecting  my 
professional  obligations  if  I did  not  offer  my  personal  congratulations 
for  a job  excellently  done  to  Mr.  E.  J.  Mock,  who  served  as  the  senior 
engineer  on  this  project.  Mr.  Mock  almost  single-handedly  developed 
the  PAREV  program  and  worked  diligently  to  ensure  continuum  among  the 
various  phases  of  this  project. 

I should  also  like  to  thank  Mr.  G.  I.  Haas,  who  supported 
Mr.  Mock  during  the  developmental  phase  of  EQSOLVE.  Last,  but  not 
least,  I thank  Mr.  R.  K.  Ellis,  the  senior  research  engineer  in  my 
section,  for  bringing  much  needed  experimental  testing  expertise  to 
the  program. 


